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ABSTRACT 



We present a new data set of transit observations of the TrES-2b exoplanet taken in spring 2009, using the 1.2m Oskar-Liihning 
telescope (OLT) of Hamburg Observatory and the 2.2m telescope at Calar Alto Observatory using BUSCA (Bonn University 
Simultaneous CAmera). Both the new OLT data, taken with the same instrumental setup as our data taken in 2008, as well as the 
simultaneously recorded multicolor BUSCA data confirm the low inclination values reported previously, and in fact suggest that the 
TrES-2b exoplanet has already passed the first inclination threshold (!„,„,! = 83.417°) and is not eclipsing the full stellar surface any 
longer. Using the multi-band BUSCA data we demonstrate that the multicolor light curves can be consistently fitted with a given set of 
limb darkening coefficients without the need to adjust these coefficients, and further, we can demonstrate that wavelength dependent 
stellar radius changes must be small as expected from theory. Our new observations provide further evidence for a change of the orbit 
inclination of the transiting extrasolar planet TrES-2b reported previously. We examine in detail possible causes for this inclination 
change and argue that the observed change should be interpreted as nodal regression. While the assumption of an oblate host star 
requires an unreasonably large second harmonic coefficient, the existence of a third body in the form of an additional planet would 
provide a very natural explanation for the observed secular orbit change. Given the lack of clearly observed short-term variations 
of transit timing and our observed secular nodal regression rate, we predict a period between approximately 50 and 100 days for a 
putative perturbing planet of Jovian mass. Such an object should be detectable with present-day radial velocity (RV) techniques, but 
would escape detection through transit timing variations. 

Key words. Stars : planetary systems - Techniques : photometry 



1. Introduction 



As of now more than 400 extrasolar planets have been detected 
around solar- like stars. In quite a few cases several planets have 
been detected to orbit a given star, demonstrating the existence of 
extrasolar planet systems in analogy to our solar system. Just as 
the planets in our solar system interact gravitationally, the same 
must apply to extrasolar planet systems. Gravitational interac- 
tions are important for the understanding of the long-term dy- 
namical stability of planetary systems. The solar system has been 
around for more than four billion years, and the understanding of 
its dynam ical stability over that period of time is still a challenge 
ISimon et al.l (,1994). In analogy, extrasolar planet systems must 
be dynamically stable over similarly long time scales, and most 
stability studies of extrasolar planet systems have been directed 
towards an understanding of exactly those long time scales. Less 
attention has been paid to secular and short-term perturbations 
of the orbit, since such effects are qu ite difficult to detect ob- 
servationally. iMiralda-Escudel (l2002b gives a detailed discus- 
sion on what secular effects might be derivable from extrasolar 
planet transits; in spectroscopic binaries the orbit inclination / 
can only be deduced in conjunction with the companion mass, 
and therefore the detection of an orbit inclination change is vir- 
tually impossible. Short-term Transit timing variations (TTVs) 
have been studied by a n umber of authors (iHolman & Murravl 
( l2005h : lAgol et al.l (l2005h ). however, a detection of such effects 
has remained elusive so far. In principle, a detection of orbit 
change would be extremely interesting since it would open up 
entirely new diagnostic possibilities of the masses and orbit ge- 



ometries of these systems; also, in analogy to the discovery of 
Neptune, new planets could in fact be indirectly detected. 

The lightcurve of a transiting extrasolar planet with known 
period allows very accurate determinations of the radii of star 
and planet (relative to the semi-major axis) and the inclination of 
the orbital plane of the planet with respect to the line of sight to- 
wards the observer. Clearly, one does not expect the sizes of host 
star and planet to vary on short time scales, however, the pres- 
ence of a third body can change the orientation of the orbit plane 
and, hence, lead to a change in the observed inclination with re- 
spect to the celestial plane. The TrES-2 exoplanet is particularly 
interesting in this context. It orbits its host star once in 2.47 days, 
which itself is very much solar-like with parameters consistent 
with solar parameters; its age is considerable and, correspond- 
ingly, the star rotates quite slowly. Its close-in planet with a size 
of 1 .2 Rjiip is among the most m a ssive known transitin g extraso- 
lar planets (IHolmanitiD (EQQ3);[Soiii^^ 

What makes the TrES-2 exoplanet orbits even more inter- 
esting, is the fact that an apparent inclination change has been 
report ed by some of us in a previous paper (iMisIis & SchmittI 
(|2009|) . henceforth called Paper I). The authors carefully mea- 
sured several transits observed in 2006 and 2008. Assuming a 
circular orbit with constant period P, the duration of an extraso- 
lar planet transit in front of its host star depends only on the stel- 
lar and planetary radii, Rs and Rj,, and on the inclination / of the 
orbit plane w.rt. the sky plane. A linear best fit to the currently 
available inclination measurements yields an apparent inclina- 
tio n decrease of 5.1 x I 0~^ deg/day. The transit mod eling both 
bv lHoIman etal.1 (l2007h and lMisIis & SchmittI (l2009f) shows the 
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Table 1. Observation summary 



Date 


Instrument 


Filter Airmass 




S eeing 


11/04/2009 


OUT 


I 1.7877- 


1.0712 


2.94" 


28/05/2009 


BUSCA 


v,b,y,l 1.8827- 


1.0453 


3.09" 



transit of TrES-2b in front of its host star to be "grazing". In 
fact, according to the latest modeling the planet occults only a 
portion of the host star and transits are expected to disappear in 
the time frame 2020 - 2025, if the observed linear trend contin- 
ues. This "grazing" viewing geometry is particularly suitable for 
the detection of orbital changes, since relatively small changes 
in apparent inclination translate into relatively large changes in 
eclipse dura t ion. A t the same time, a search f or possible TTVs by 
iRabus et al] (|2009|) has been negative; while iRabus et all (120091) 
derive a period wobble of 57 sec for TrES-2, the statistical qual- 
ity of their data is such that no unique periods for TTVs can be 
identified. 

The purpose of this paper is to present new transit observa- 
tions of the TrES-2 exoplanet system obtained in 2009, which 
are described in the first part (Sec|2]i and analysed (Sec. [Sj. In 
the second part (Sec|4]i of our paper we present a quantitatively 
analysis of what kind of gravitational effects can be responsi- 
ble for the observed orbit changes of TrES-2b and are consistent 
with all observational data of the TrES-2 system. 

2. Observations and data reduction 

We obser ved two transits of TrES -2 usin g the ephemer i s sug- 
gested bv lO'Donovan et all (l2006l) and by iHolman eTal] (120071) 
from = 2, 453, 957.6358 + E ■ (2.47063 days), 

using the 1.2m Oskar Liihning telescope (OLT) at Hamburg 
Observatory and Calar Alto Observatory 2.2m telescope with 
BUSCA. 

The OLT data were taken on 11 April 2009 using a liKxiK 
CCD with a lO'xlO' FOV and an I-band filter as in our previous 
observing run (Paper I). The readout noise and the gain were 
16.37 and 1.33 /ADU respectively. With the OLT we used 
60-second exposures which provided an effective time resolu- 
tion of 1.17 minutes. During the observation, the airmass value 
ranged from 1.7877 to 1.0712 and the seeing was typically 2.94". 

The Calar Alto data were taken on 28 May 2009 using 
BUSCA and the 2.2m telescope. BUSCA is designed to per- 
form simultaneous observations in four individual bands with 
a FOV of ll'xir. Therefore it has four individual AKxAK 
CCD systems which cover the ultra-violet, the blue-green, the 
yellow-red and the near-infrared part of the spectrum (chan- 
nel a-d respectively). For this run we used the Stromgren fil- 
ters V (chn. a), b (chn. b), and y (chn. c) , and a Cousins-I fil- 
ter for the near-infrared (chn. d). The readout-noise for these 
four detectors are 9.09e", 3.50e", 3.72e", and 3.86e" respec- 
tively for the a,b,c, and d channels. The gain values for the same 
channels are \3Ale'IADU, l.lble-jADU, I.133e-/ADU, and 
1.73le^ /ADU respectively. The airmass value ranged from 
1.8827 to 1.0453 and the seeing was 3.09". For the BUSCA 
observations we took 30 seconds exposure yielding an effective 
time resolution of 1.63 minutes. In Table [T] we summarize the 
relevant observation details. 

For the data reduction, we used Starlink and DaoPhot soft- 
ware routines, and the MATCH code. We perform the normal re- 
duction tasks, bias subtraction, dark correction, and flat fielding 
on the individual data sets before applying aperture photometry 



Table 2. Relative Photometry data. 



HJD Relative Flux Uncertainty Flag 



2454933.44031 


0.99172 


0.0037 


OLT - 


- 1 


2454933.44091 


0.99204 


0.0038 


OLT - 


- 1 


2454933.44221 


1.00514 


0.0039 


OLT - 


- 1 


2454933.44281 


1.00314 


0.0038 


OLT - 


-1 


2454933.44341 


0.99458 


0.0035 


OLT - 


- 1 


2454933.44401 


0.99745 


0.0037 


OLT - 


-I 



Note : Table 2: Relative photometry vs. time; note that the complete 
table is available only electronically. The time stamps refer to the 
Heliocentric Julian Date (HJD) at the middle of each exposure. The 
'Flag' column refers to the telescope and filter. 

on on all images. For TrES-2, we selected the aperture photome- 
try mode using apertures centered on the target star, check stars, 
and sky background. Typical sky brightness values for the 1 1 
April and for the 28 May were 89 and 98 ADUs, respectively, 
i.e., values at a level 0.008% and 0.006% of the star's flux, re- 
spectively (for I-filter). For the relative photometry we used the 
star U 1350- 10220 179 as a reference star to test and calibrate the 
light curve. For the data analysis presented in this paper we did 
not use additional check stars, but note that we already checked 
this star for constancy i n the Paper I. To estima te the magni- 
tude errors, we followed lHowell & Everetl (1200 ll) and the same 
procedure as described in our first paper (Paper I) to obtain bet- 
ter relative results. Our final relative photometry is presented in 
Table |2] which is available in its entirety in machine-readable 
form in the electronic version of this paper. 

3. IVIodel analysis and results 

In our model analysis we proceeded in exactly the same fashion 
as described in Paper I. Note that the assumpt ion of circularit y 
appears to apply very we ll to TrES-2b (.O' Donovan et al.l (l2009l) ; 
iO' Donovan et al.1 (l2006h ); the assumption of constant period and 
hence constant semi-major axis will be adressed in section 4. For 
our modelling we specifically assumed the values - 1 .0037?o, 
Rp = 1.222 Rj, P = 2.470621 days, a = 0.0367 AU for the radii 
of star and planet, their period and the orbit radii is respectively . 
All limb darkening coefficients were taken from IClare3 ( |2004|) . 
and for the OLT data we used the same values as in Paper I, viz.. 
Mi =0.318 and U2 - 0.295 for the I filter, as denoted by SI in 
Tabled 

3.1. OLT data & modeling 

Our new OLT data from 1 1 April 2009 were taken with the same 
instrumental setup as our previous data taken on 18 September 
2008. The final transit light curve with the best fit model is 
shown in Fig[T] the reduced light curve data are provided in 
Table [1] Keeping the planetary and stellar radii and the limb 
darkening coefficients fixed (at the above values), we determine 
the best fit inclination and the central transit time Tc with our 
transit model using the x' method; the thus obtained fit results 
are listed in Table |3] The errors in the derived fit parameters are 
assessed with a bootstrap method explained in detail in Paper 
I, however, we do not use random residuals for the new model, 
but we circularly shift the residuals after the model substraction 
to produce new light curves for the bootstrapped data following 
Alonso et al. (2008). In Tab. [3] we also provide corrected cen- 
tral times and errors for those transists already reported in pa- 
per I, since due to some typos and mistakes the numbers quoted 
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Fig. 1. Observed TrES-2 light curves and model fits for the light 
curve taken with the 1 .2m OLT at Hamburg Observatory taken 
in April 2009. 
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Fig. 2. Four light curves and corresponding model fits taken in 
May 2009 with the 2.2m telescope at Calar Alto using BUSCA. 



for central time and their error are unreliable (five first lines of 
Table [3]). In Fig. [3] we plot the thus obtained inclination angle 
distribution for epoch April 2009 (solid curve) compared to that 
obtained in September 2008 (dashed curve). While the best fit in- 
clination of / = 83.38° differs from that measured in September 
2008 (i = 83.42°), the errors are clearly so large that we cannot 
claim a reduction in inclination from our OLT data (taken with 
the same instrumental setup, i.e, in September 2008 and April 
2009) alone, however, is clearly consistent with such a reduc- 
tion. 




Fig. 3. Inclination distribution for the OLT data using 1000 
Monte Carlo simulations for the September 2008 (solid curve) 
and April 2009 (dashed curve); see the text for details. 

Table 3. Individual values of duration, inclination, values and 
limb darkening coefficients from five light curve fits; units for 
duration and errors are minutes and for inclination and errors 
are degrees. The OLT light curve has 268 points, the I-filter light 
curve has 198, the y-filter light curve has 208, the b-filter light 
curve 198 and the v-filter Ught curve 162 points respectively. 



Tc time [HJD] Duration Errors Inclin. Errors X'^value LDL 

3989.7529 ± 0.00069 110.308 0.432 83.59 0.019 4321 

3994.6939 ± 0.00066 109.230 0.448 83.56 0.019 296.8 

4041.6358 ±0.00070 109.025 0.430 83.55 0.019 290.6 

4607.4036 ± 0.00072 106.620 0.883 83.44 0.036 179.1 

4728.4640 ±0.00071 106.112 0.870 83.43 0.036 190.1 

4933.5274 ± 0.00076 105.077 0.964 83.38 0.039 296.1 51 

4980.4675 ± 0.00068 105.056 0.848 83.38 0.034 181.7 51 

4980.4679 ± 0.00090 104.748 1.076 83.37 0.043 202.2 52 

4980.4667 ± 0.00082 103.832 1.021 83.33 0.041 186.0 53 

4980.4678 ± 0.00100 104.363 1.194 83.35 0.048 195.0 54 

4980.4675 ± 0.00060 104.522 83.36 0.030 



3.2. BUSCA data and modeling 

Our BUSCA data have the great advantage of providing 
simultaneously measured multicolor data, which allows us 
to demonstrate that limb darkening is correctly modelled and 
does not affect the fitted inclinations and stellar radii. The limb 
darkening coefficients used f or the analysis of the BUSCA data 
were also taken from IClaretl (|2004|) ; we specifically used ui - 
0.318 and U2 - 0.295 for the I filter (SI), i.e., the same values as 
for the OLT, ui = 0.4501 and U2 = 0.286 for the y filter (S2), ui 
= 0.5758 and U2 = 0.235 for the b filter (S3), mi = 0.7371 and 
M2 = 0. 1117 for the v filter (S4) respectively. The data reduction 
and analysis was performed in exactly the same fashion as for 
the OLT data, we also used the same comparison star U1350- 
10220179; the reduced light curve data are also provided in 
Table |2] The modelling of multicolor data for extrasolar planet 
transits needs some explanation. In our modelling the host star's 
radius is assumed to be wavelength independent. Since stars do 
not have solid surfaces, the question arises how much the stellar 
radius does actually change with the wavelength. This issue 
has been extensively studied in the solar context, where the limb 
of the Sun can be directly observed as a function of wavelength 
(lNeckellll995h . Basically the photospheric height at an optical 
depth of unity is determined by the ratio between pressure 
and the absorption coefficient, and for the Sun, Neckel (19951) 
derives a maximal radius change of 0.12 arcsec between 3000 
and 10000 A, which corresponds to about 100 km. If we assume 
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Table 4. tests for two different inclination values and the 
errors after 1000 Monte Carlo simulations. 



Inclination [°] 


OLT 


B-I 


B-y 


B-b 


B-v 


83.57 (Holman) 


304.9 


199.4 


207.8 


216.0 


206.3 


cr^i (Holman) 


16.12 


14.04 


16.16 


21.05 


14.26 


83.36 {this paper) 


213.8 


173.44 


164.3 


180.0 


131.1 


(T^i {this paper) 


11.96 


11.11 


9.96 


12.01 


8.77 



similar photospheric parameters in the TrES-2 host star, which 
appears reasonable since TrES-2 is a GOV star, we deduce a 
wavelength-dependent radius change on the order of 100 km, 
which is far below a percent of the planetary radius and thus not 
detectable. Therefore in our multi-band data modelling we can 
safely fix the radius of the star and the radius of the planet in a 
wavelength-independent fashion. 

We thus kept the stellar and planetary fixed at the above values, 
treated all BUSCA channels as independent observations and 
fitted the light curves - as usual - by adjusting inclination 
and central transit time. The filter light curves are shown 
together with the so obtained model fits in Fig. |2l the model 
fit parameters are again summarized in Table [3] We emphasize 
that we obtain good and consistent fits for all light curves 
with the chosen set of limb darkening coefficients, thus demon- 
strating our capability to correctly model multicolor light curves. 

Since the BUSCA data are recorded simultaneously, it is 
clear that the light curves in the four BUSCA channels must ac- 
tually be described by the same values of inclination and tran- 
sit duration. We therefore simultaneously fit all the four light 
curves, leaving free as fit parameters only the central time 
and inclination /. With this approach we find an average inclina- 
tion of / - 83.36° ± 0.03°, which is consistent with our spring 
OLT data and also suggests that the inclination in spring 2009 
has further decreased as compared to our 2008 data. 



0.4 7 




83.20 53.25 83.30 83.35 83.40 83.45 33.50 

hclination ["] 

Fig. 4. Inclination distribution derived 1000 Monte Carlo sim- 
ulations for the multi-band BUSCA observations data in the 
I (crosses), y (stars), b (diamonds) and v (x-symbol) filter 
bands; the mean inclination distribution from the four BUSCA 
lightcurves is shown as solid line, the inclination distribution de- 
rived from the OLT data taken in April 2009 is also shown for 
comparison (dashed line). 



3.3. Joint modeling 

Using all our data we can check whether our assumption of the 
wavelength-independence of the radius is consistent with the ob- 
servations. For this consistency test we kept the inclination value 
fixed at / = 83.36° and fitted only the radius of the star Rs and 
Tc. The errors on Rg were again assessed by a Monte Carlo sim- 
ulation as described in Paper I and the distribution of the thus 
derived stellar radius values Rs is shown in Fig. 21 as apparent 
from Fig. |4] all BUSCA channels are consistent with the same 
stellar radius as - of course - expected from theory since any 
pulsations of a main-sequence star are not expected to lead to 
any observable radius changes. 

The crucial issue about the TrES-2 exoplanet is of course the 
constancy or variability of its orbit inclination. Our new BUSCA 
and OLT data clearly support a further decrease in orbit inclina- 
tion and hence decrease in transit duration. In order to demon- 
strate the magnitude of the effect, we performed one more se- 
quence of fits, this time keeping all physical parameters fixed 
and fitting only the central transit times Tc using 1000 Monte 
Carlo realisations and studying the resulting distribution in;^'^; 
for the inclination we assumed for, first, the value / = 83.57° 
(as derived by Holman et al. for their 2006 data) and the value 
/ = 83.36° (this paper for the 2009 data). The fit results (in terms 
of obtained values) are summarised ed in Tab.|4] which shows 
that for all (independent) data sets the lower inclination values 
yield smaller ;t'^-values; for some filter pairs the thus obtained 
improvement is extremely significant. 

0-25 




Star Radius [R./RJ 

Fig. 5. Stellar Radius distribution (/?»), derived from 1000 Monte 
Carlo simulations in four different filters (from higher to lower 
peak, I, y, b and v filter respectively). The overlap of the curves 
shows that all colors can be ex plained with th e same stellar ra- 
dius as suggested by theory (see iNeckell d 1 995h ) . 



3.4. Inclination changes 

In Fig. |6] we plot all our current inclination vs. epoch data to- 
gether with a linear fit to all data; a formal regression analy- 
sis yields for the time evolution of the inclination (/ = ig + a ■ 
(Epoch), io = 83.5778, a = 0.00051). In Paper I we noted the 
inclination decrease and predicted inclination values below the 
first transit threshold (/,„,„_i < 83.417°) after October 2008. Both 
the new OLT data set and all BUSCA channel observations yield 
inclinations below the first transit threshold. While the error in a 
given transit light curve is typically on the order of 0.04° for /, 
we consider it quite unlikely that 5 independent measurements 
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all yield only downward excursions. We therefore consider the 
decrease in inclination between fall 2008 and spring 2009 as 
significant, conclude that the inclination in the TrES-2 system 
is very likely below the first transit threshold, and predict the 
inclinations to decrease further; also, the transit depths should 
become more and more shallow since the exoplanet eclipses less 
and less stellar surface. 




Epoch [E] 

Fig. 6. Epoch versus inclination together with a linear fit to the 
currently available data ; the diamond points are those taken in 
2006 by Holman et al. (2007), and those taken in 2008 and re- 
ported in Paper 1. The square points are derived from our new 
observations taken in April and May 2009. The solid lines show- 
ing two linear fits, from the first paper (dashed line) and the fit 
from the present paper (solid line). 



3.5. Period changes 



The observed change in orbit inclination is in marked contrast 
to the period of TrES-2b. While po ssible TTVs in TrES-2b have 
been studied bv lRabus et al.l(l2009l) we investigate the long-term 
stability of the period of TrES-2b. From our s even tran s it mea - 
surements (plus five more data points of Rab us et all ( |2009l) ) 
spanning about ~ 400 eclipses we created a new O-C diagram 
(cf.. Fig. [7]i; note that w e refra ined from using the transit times 
dicussed bv lRaetzetaT] ( l2009h . since these transits were taken 
with rather inhomogeneous instruments and sometimes only par- 
tial transit coverage. For our fit we used a modified epoch equa- 
tion HJDc = HJDo + E P, where we set P = Po+ Pit - HJDo) 
and explicitly allow a non-constant period P. We apply a fit 
to find the best fit values for P, P,,, and HJDo- With this ap- 
proach we find a best fit period change of P = 5 x 10"'^, how- 
ever, carrying out the same analysis keeping a fixed period shows 
that the fit improvement due to the introduction of a non-zero P 
is insignificant. We then find as best fit values Po - 2.47061 
and HJDq — 2453 957.6350 conforming to the values derived by 
iRabus et"aLl (12009). Thus, the period of TrES-2b is constant, and 
any possible period change over the last three years must be less 
than about 1 sec. 



4. Theoretical implications of observed inclination 
change 

The results of our data presented in the preceding sections 
strenghten our confidence that the observed inclination changes 
do in fact correspond to a real physical phenomenon. Assuming 



200 
Epo 



Fig. 7. O-C values versus epoch including the transits observed 
by Holman et al., Rabus et al. and our data denoted by triangles, 
squares and diamonds respectively. 



now the reality of the observed inclination change of A/ ~ 
0.075°/yr, given the constancy of the period and the absence 
of TTVs at a level of a; 100 sec we discuss in the following a 
physical scenario consistent with these observational findings. 
We specifically argue that the apparent inclination change should 
be interpreted as a nodal regression and then proceed to examine 
an oblate host star and the existence of an additional perturb- 
ing object in the system as possible causes for the change of the 
orbital parameters of Tres-2b. 



4. 1. Inclination change or nodal regression? 

It is important to realize that the reported apparent change of 
inclination refers to the orientation of the TrES-2 orbit plane 
with respect to the observer's tangential plane. It is well known 
that the z-component of angular momentum for orbits in an az- 
imuthally symmetric potential is constant, resulting in a constant 
value of inclination. An oblate star (cf.. Sec. 14.2b or the aver- 
aged potential of a third body (cf.. Sec. 14.31 ) naturally lead to 
such potentials, thus yielding orbits pre cessing at (more or less ) 
constant inclination as also realized bv lMiralda-Escudel ( l2002h . 
Such a precession would cause an apparent inclination change, 
however, physically this would have to be interpreted as nodal 
regression at fixed orbit inclination. To interpret the observa- 
tions, one has to relate the rate of nodal regression to the rate 
of apparent inclination change. 

Consider a (massless) planet of radius Rp orbiting a star of 
radius at some distance d. Let the planet's orbit lie in a plane 
with a fixed inclination / relative to the x-y plane, which we take 
as invariant plane. Let an observer be located in the x-z plane 
with some elevation y, reckoned from the positive x-axis. The 
line of the ascending node in the x-y plane is denoted by the 
angle Q., with = implying the ascending node pointed along 
the negative y-axis (see Fig. [8). Let in the thus defined geometry 
the angle denote the angle between planet and observer as seen 
from the central star. For each system configuration defined by 
the angles (y, /, Q) there is a minimal angle ^'min between orbit 
normal and observer obtained in each planetary orbit, which can 
be computed from 

COS^Ftnin = «Orbit ' ''Gbserver - - COS y sin / COS Q -I- sinyCOS/. (1) 



A transit takes place when 
\cos'¥,^J <(Rp+R,)/d, 



(2) 
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Fig. 8. System geometry (see text for details): The orbit plane 
with normal vector noibit is inclined relative to the x-y plane by 
an inclination the observer (towai^ds robserver) is within the x-z 
plane. 



and from the geometry it is clear that the observed inclination 
Jobs, i-c, the parameter that can be derived from a transit light 
curve is identical to *Pmin- Setting then ^F^in = 'obs and differen- 
tiating Eq. [T]with respect to time we obtain 



di'obs sin n cos y sin / dQ 



df 



sin iobs 



df 



(3) 



Eq. [3] relates the nodal regression of the orbit to its corre- 
sponding observed apparent rate of inclination change d/obs/df, 
given the fixed inclination / relative to the x-y plane. Since tran- 
sit observations yield very precise values of /obs, the required as- 
cending node Q. and its rate of change can be computed, once the 
orbit geometry through the angles y and /obs is specified. In Fig.|9] 
we show a contour plot of the linear coefficient of Eq.[3]between 
nodal regression and observed inclination change as a function 
of orbit geometry. Note that the apparent change of inclination 
due to the nodal regression does vanish for / = 0. Physically it 
is clear that a perturbing planet in a coplanar orbit cannot exer- 
cise a torque and therefore cannot cause the observed inclination 
variation. We will therefore always assume ; '^^^ in the follow- 
ing. 

4.2. Ob/afe host star 

We first consider the possibility that the TrES-2 host star is 
oblate. The motion of a planet around an oblate host star is equiv- 
alent to that of an artifical satellite orbiting the Earth, a prob- 
lem intensely s t udied over the last decades and well understood 
IConnon Smith (2005"). The potential U(r, (p) of an axisymmetric 
body of mass M and radius R can be expressed as a power series 
involving the so-called harmonic coefficients. In second order 
one approximates the potential U (r, (p) as 



f/(r,0) 



GM 



1 



3 sin" - l) 



(4) 



where r is the radial distance from the body's center, is lati- 
tude above the equator and G denotes the gravitational constant. 
Clearly, the perturbing term in the potential is proportional to J2 
and a perturbation calculation yields as first order secular pertur- 
bation the angular velocity of the ascending node as 



dQ 
d7 



3 JiR^ 



2(1 -e2)2 



cos I 



7/3 



1 



(GM,otal)2/3' 



(5) 
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Fig. 9. Linear coefficient of Eq. [3]between nodal regression and 
observed inclination change as a function of view geometry (cf., 
Eq. [3j, computed for /obs = 83.38° as applicable for TrES-2; 
the plotted contour levels denote values of 1.1, 1.3, 1.5, 2., 2.5, 
5. ,4. ,6. ,8. from right bottom up. 



where e and / denote the eccentricity and inclination of the orbit- 
ing body, P its period, Mtotai the sum of the masses of planet and 
host star and the validity of Kepler's third law has been assumed. 

Interpreting the observed inclination change as nodal re- 
gression due to an oblate host star, we can compute a lower 
bound on the required harmonic coefficient J2 assuming e - 0. 
Therefore, we combine Eqs.[3]and|5]to obtain an expression for 
J2- Excluding pathological cases like i,D. = and j - njl and 
neglecting the planetary mass in Eq.|5]we find for any given set 
of parameters (Q, y, i) 



Ji > JT = 



2 di'obs 
3~dr 



-1 
27r/ 



7/3 



(GMhost)'^' 



R^ 



■ sin iobs, 



(6) 



where Mhost denotes the mass of the host star. Mass, radius, incli- 
nation and period of the TrES-2 system are well known, and us- 
ing the measured nodal regression we find 72,TiES-2 ~ 1 .4 x 10""*, 
i.e., a value smaller than that of the Earth (72,® - 0.00108) by an 
order of magnitude, but considerably large r than that of the Sun , 
which is usually taken as 72,o ~ 3-6x10'^ iRozelot et al J ( 1200 lb . 
Since the host star of TrES-2 is a slow rotator very similar to the 
Sun Sozzetti et al. (2007), we expect similarly small J2 values in 
contrast to our requirements. We therefore conclude that oblate- 
ness of the host star cannot be the cause for the observed orbit 
variations. 

4.3. Perturbation by a third body 

An alternative possibility to explain the observed orbit variations 
of the TrES-2 exoplanet would be the interaction with other plan- 
ets in the system. Let us therefore assume the existence of such 
an additional perturbing planet of mass nip, circling its host star 
of mass mo with period Pp at distance rp located further out com- 
pared to the known transiting TrES-2 exoplanet. This three-body 
probl em has been con s idered in p ast in the context of triple sys- 
tems dKhafiullin et all d 199 lb . 10 (120061) ) and the problem of ar- 
tificial Earth satellites, whose orbits are perturbed by the Moon. 
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In lowest order, the perturbing gravitational potential R2 onto the 
inner planet with mass m and distance r is given by the expres- 



mp + m„\Pp ) \rpj \al 



r \2 3 cos^ 5-1 



(7) 



where the angle S denotes the elongation between the perturbed 
and perturbing planet as seen from the host star, a and Up de- 
note their respective semi-major axes and the validity of Kepler's 
third law has been assumed. Note that in this approach the per- 
turbed body is assumed to be massless, implying that its pertur- 
bations onto the perturbing body are ignored. Next one needs to 
insert the orbital elements of the two bodies and, since we are 
interested only in secular variations, average over both the peri- 
ods of the perturbed and p erturbing planet . This is the so-called 
double-averaging method jBrouckel d2003l) ). which, howe ver, in 
mor e or less the same form has also been applied by Q ( l2006l) 
and iKovalevskvl (|1967|) . Denoting by e the eccentricity of the 
perturbed planet, by lj the longitude of the periastron and by / 
the angle between the tw o orbital planes, one obtains after some 
lengthy computation (see IKovalevskvl (1 1 9671) ) 

m„ I na \^ 

Ri^ —\TB-] xKo{i,e,aj), (8) 

nip + nio \2Ppj 

with the auxiliary function KQ{i, e, ai) given by 

Ko{e, i, oj) - (6 cos^ / - 2) -1- ^9 cos" / - 3) -1- I5e^ sin^ / cos 2lo. 

The partial derivatives of R2 with respect to the orbital elements 
are needed in the so-called Lagrangian planetary equations to de- 
rive the variations of the orbital elements of the perturbed body. 
One specifically finds for the motion of the ascending node 



dQ nip 3n P cos/ 



X Ki{e,at), 



dt nip + nio 4 P^ y/i _ g2 
where the auxiliary function Ki (e, cS) is defined through 
K\{e, (jj) - 5e^ cos 2ci) - 3e' - 2 
and for the rate of change of inclination 



di 2 '^p 15;r P 1 
dt nip + nio 8 Pi VT^^e^ 



sin (2/) sin(2(jj). 



(9) 



(10) 



(11) 



As is obvious from Eq. [TTj the rate of change of inclination in 
low eccentricity systems is very small and we therefore set e -0. 
Assuming next a near coplanar geometry, i.e., setting cos / x 1, 
we can simplify Eq.|9]as 



dO. nip 37T P 

dt nip + m„ 2 P\ 



(12) 



If we assume a host star mass and interpret the observed incli- 
nation change as the rate of nodal regression via Eq. [3] Eq. [T2| 
relates the unknown planet mass nip to its orbital period Pp. 

4.3.1 . Sanity check: Application to tine Solar System 

The use of Eq. [12] involves several simplifications. Thus, it is 
legitimate to ask, if we are justified in expecting Eq. [T2]to de- 
scribe reality. As a sanity check we apply Eq. [12] to our so- 
lar system. Consider first the motion of the Moon around the 
Earth, i.e., P - 27.3 d, which is perturbed by the Sun, i.e.. 



Pp - 365.25 d. Since for that case our nomenclature requires 
nip » Too and / ~ 5.1°, we find from Eq. [12] a time of 17.83 
years for the nodes to complete a full circle, which agrees well 
with the canonical value of 18.6 years for the lunar orbit. In the 
lunar case it is clear that the Sun with its large mass and close 
proximity (compared to Jupiter) is by far the largest perturber 
of the Earth-Moon two-body system and this situation is exactly 
the situation described by theory. 

Consider next the the perturbations caused by the outer 
planets of our solar system. Considering, for example, Venus, 
we can compute the perturbations caused by the planets Earth, 
Mars, Jupiter, Saturn and Uranus. Since the perturbation strength 
scales by the ratio nipPp^, we can set this value to unity for the 
Earth and compute values of 0.03, 2.26, 0.1 1 and 0.002 for Mars 
through Uranus respectively. So clearly, Venus is perturbed by 
several planets, but the perturbations by Jupiter are strongest. We 
therefore expect that our simple approach is not appropriate. We 
further note that among the outer solar system planets long pe- 
riod perturbations and resonances occur, which are not described 
by Eq. [12] If we nevertheless compute the nodal regression for 
Venus caused by Jupiter using Eq.[T2] we find a nodal regression 
of 0.1°/cty for Venus, and 0.3°/cty for Mars. Using the orbital 
elements computed by ISimon et al.l ([1994) and calculating the 
nodal regressions of Venus and Mars in the orbit plane of Jupiter 
we find values smaller than the true values, but at least, they 
computed values in the right order of magnitude and do not lead 
to an overprediction of the expected effects. 



4.3.2. Sanity check: Application to V 907 Sco 



From archival studies iLacv et all d 19991) report the existence of 
transient eclipses in the triple star V 907 Sco. According to 
iLacv et all (1 19991) this system is composed out of a short-period 
(Pshort = 3.78 days) binary containing two main sequence stars 
of spectral type ~ AO and mass ratio 0.9, orbited by a lower mass 
third companion (Piong = 99.3 days), of spectral type mid-K or 
possibly even a white dwarf. The close binary system showed 
eclipses from the earliest reported observations in 1899 until 
about 1918, when the eclipses stopped; eclipses reappeared in 
1963 and were observed until about 1986. I nterpreting t he ap- 
pearance of eclipses due to nodal regression, ILacv etal.1 (119991) 
derive a nodal period of 68 years for V 907 Sco. Using Eq. [12] 
and assuming a mass of 2 Mq for the host (m„) and a mass of 
0.5 Mq for the companion (nip), we compute a nodal regression 
period of 4 7.6 years, which a grees well with the nodal period es- 
timated bv lLacv et al.l ( Il999l) . We therefore conclude that Eq.[T2] 
also provides a reasonable description of the transient eclipse 
observations of V 907 Sco. 



4.3.3. Application to TrES 2 

Applying now Eq. [12] to the TrES-2 exoplanet we express the 
period of the (unknown) perturbing planet as a function of its 
also unknown mass through 



37rP 1 



2 dQ./dt 



(13) 



Since the host star mass and the nodal regression are known, 
the perturbing mass is the only remaining unknown; we note that 
Eq. [13] should be correct, as long as there is only one domi- 
nant perturber in the system with a low-eccentricity orbit suf- 
ficiently far away from the known close-in transiting planet. In 
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Fig. [To] we plot the expected perturber period Pp as a function 
of lUp (in Jupiter masses) for m,, - 1 and the measured 
rate of nodal regression assuming a linear coefficient of unity 
in Eq. [3] Assuming ad hoc a mass of about one Jupiter mass for 
this perturber and taking into account that a factor of a few is 
likely (cf., Fig.|9l), we find that periods of 50 to 100 days are re- 
quired to explain the observed nodal regression in TrES-2. Such 
an addition al planet should be rel atively easily detectable with 
RV studies. [Paemgen et alj (l2009h report the presence of a faint 
companion about one arcsec away from TrES-2. Assuming this 
companion to be physical, a spectral type between K4.5 and K6, 
a mass of about 0.67 at a distance of 230 AU with a pe- 
riod of 3900 years follow. Computing with these numbers the 
maximally expected nodal regression from Eq.[T2] one finds val- 
ues a couple of orders of magnitude below the observed values. 
We thus conclude that this object cannot be the cause of the ob- 
served orbit variations. On the other hand, this companion, again 
if physical, makes TrES-2 particularly interesting because it pro- 
vides another cases of a planet/planetary system in a binary sys- 
tem, and eventually the orbit planes of binary system and the 
planet(s) can be derived. 
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Fig. 10. Period of hypothesized second planet vs. mass assuming 
a linear coefficient of unity in Eq.[3]and near coplanarity. 



4.4. Transit timing variations by a putative perturber 

A perturbing second planet capable of causing fast nodal pre- 
cession on the transiting planet is also expected to induce short- 
term periodic variations of its orbital elements. In addition to the 
secular precession of the node of the orbit we would thus ex- 
pect to see short-term transit timing variations (TTVs), periodi c 
variations of the mid-transit times (iHolman & Murravi (l2005h : 
lAeol et al.l(l2005h '). Just as nodal regression, the TTV signal can 
be used to find an d characterize plan etary companions of tran- 
siting exoplanets. iRabus et"aLl (l2009l) carefully analyzed eight 
transit light curves of TrES-2 over several years. However, they 
were unable to detect any statistically significant TTV ampli- 
tudes in the TrES-2b light curves above about 50 s; cf. their 
Fig. 10. Therefore, the existence of perturbing objects leading 
to TTVs on the scale of up to about 50 s are consistent with 
actual observations. Putting it diff'erently, the orbital parameters 
of any perturbing object causing the nodal precession of the or- 
bit should yield a TTV amplitude below that and hence remain 
undetactable in the presently existing data. 

To analyze the mutual gravitational influence of a perturb- 
ing second planet in the system on TrES-2, we have to treat the 



classical three-body problem of celestial mechanics. Instead of 
direct A^-body integrations of the equations of motion we use an 
alternative method based on anal ytic perturbation theory devel- 
ope d and extensivel y tested by iNesvornv & Morbidellii (l2008l) 
and lNesvornvl (|2009|) . Outside possible mean-motion resonances 
their approach allows for a fast computation of the expected TTV 
amplitude given a combination of system parameters. As input 
we have to specify the orbital elements and masses of both plan- 
ets. Consistent with the observations we assume TrES-2 to be in 
a circular orbit around its host star, while we allow for diff'er- 
ent eccentricities Cp and periods Pp of the perturber, which we 
assume to be of Jovian mass; since the TTV amplitudes scale 
nearly linearly with the perturber mass, we confine our treat- 
ment to nip = 1 Mjup; all other orbital elements are set to zero. 
This is justified as these parameters in most cases d o not lead 
to a significant amplification of the TTV signal (see iNesvorn^ 
(I2009h for a detailed discussion of the impact of these orbital el- 
ements). The resulting TTV amplitude for different reasonable 
orbit configurations (given the observed secular node regression 
as discussed above) of the system of Ppci = 30, 50 and 70 days 
is plotted vs. the assumed eccentricity in Fig. (TTJ th e currently 
available upper limit to any TTV signal derived by Rab us et akl 
(12009 ) is also shown. As is obvious from Fig.[TT] a Jovian-mass 
perturber at a distance required to impose the observed secular 
changes (period of 50 - 100 days) leads to a TTV signal well 
below the current detection limit for all eccentricities as long 
as gp < 0.4. We therefore conclude that a putative perturbing 
Jovian-mass planet with a moderate eccentricity and with a pe- 
riod between 30 - 70 days would not yield any currently de- 
tectable TTV signal and would therefore be a valid explanation 
for the observed inclination change in the TrES-2 system. 
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Fig. 11. Amplitude of expected Transit Timing Variations 
(TTVs) in the TrES-2 system. The perturber is assumed to have 
nip - I Mjup. Its eccentricity and period Pp are varied within 
plausible ranges. The orbit of TrES-2 is assumed to be circular 
The ver tical do tted line marks the best fitting TTV signal found 
bv iRabus et al.1 ((20091) of 57 s. 



5. Conclusions 

In summary, our new observations taken in the spring of 2009 
confirm the smaller transit durations reported in Paper I and sug- 
gest an even further decrease. With our simultaneously taken 
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multicolor BUSCA data we demonstrate that the recorded mul- 
ticolor lightcurves can be consistently modelled with a reason- 
able set of limb darkening coefficients, and that there is no 
need to fit the limb darkening coefficients to any particular light 
curve. An error in the description of the limb darkening there- 
fore appears thus as an unlikely cause of the observed inclination 
changes. Also as expected, the obtained stellar radius is inde- 
pendent of the wavelength band used, demonstrating the internal 
self-consistency of our modelling. 

As to the possible causes for the observed apparent orbit in- 
clination change in TrES-2 we argue that the apparent observed 
inclination change is very likely caused by nodal regression. The 
assumption of an oblate host star leads to implausibly large J2 
coefficients, we therefore favor an explanation with a third body. 
We argue that Eq. [12] is a reasonable approximation for the in- 
terpration of the observed inclination changes; applying it to the 
TrES-2 system, we find that a planet of one Jovian mass with 
periods between 50 - 100 days would suffice to cause the ob- 
served inclination changes, while at the same time yield TTVs 
with amplitudes well below the currently available upper limits. 

The assumption of such an additional planet in the TrES-2 
system is entirely plausible. First of all, if it is near coplanar 
with TrES-2b, it would not cause any eclipses and therefore re- 
main undetected in transit searches. Next, an inspection of the 
exosolar planet data base maintained at www.exoplanet.eu re- 
veals a number of exoplanet systems with properties similar to 
those postulated for TrES-2, i.e., a close-in planet together with 
a massive planet further out: In the Gl 581 system there is a 
0.02 Jupiter-mass planet with a period of 66 days, and in fact 
a couple of similarly massive planets further in with periods of 
3.1, 5.4 and 12.9 days respectively; in the system HIP 14810 
there is a close-in planet with a 6.6 day period and a somewhat 
lighter planet with a period of 147 days, in the HD 160691 sys- 
tem the close-in planet has a period of 9.6 days and two outer 
planets with Jupiter masses are known with periods of 310 and 
643 days. It is also clear that in these systems nodal regression 
changes must occur, unless these systems are exactly coplanar, 
which appears unlikely. Therefore on longer time scales the ob- 
served orbit inclination in these systems must change, but only 
in transiting systems the orbit inclination can be measured with 
sufficient accuracy. Because of its apparent inclination change 
TrES-2b is clearly among the more interesting extrasolar plan- 
ets. If the system continues its behavior in the future the tran- 
sits of TrES-2b will disappear. Fortunately, within the first data 
set of the Kepler mission ~ 30 transits should be covered. From 
our derived inclination change rate of Ai ~ 0.075°/yr this cor- 
responds to an overall change of Ai ~ 0.015° in this first data 
set, which ought to be detectable given the superior accuracy of 
the space-based Kepler photometry. As far as the detection of 
our putative second planet is concerned, RV methods appear to 
be more promising than a search for TTVs, unless the orbital 
eccentricities are very large. 
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